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Byram's  Energy  Criterion 
for  Wildland  Fires:  Units 
and  Equations 

Ralph  M.  Nelson,  Jr. 


Abstract — The  equations  presented  by  G.  M.  Byram  in  1959  for  calculating  ener- 
gy flow  rate  in  the  wind  field  (Pw)  and  in  the  convection  column  (Pf)  above  a  wild- 
land  fire  are  inconsistent  because  units  for  each  rate  are  not  those  of  an  energy  rate 
(or  power)  per  unit  area.  New  equations  for  Pw  and  Pf  are  introduced,  and  sample 
calculations  are  presented  for  the  British  gravitational  and  SI  systems  of  units. 
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The  concepts  of  energy  flow  rate  in  the  wind  field  and  in  the  convection 
column  above  a  line  of  fire  were  introduced  by  G.  M.  Byram  (1959).  The  ener- 
gy flow  rate  in  the  wind  field  (Pw)  is  a  rate  of  flow  of  kinetic  energy  through 
a  vertical  plane  of  unit  area  at  a  specified  height  in  a  neutrally  stable  atmo- 
sphere, one  in  which  the  temperature  of  unsaturated  air  decreases  9.8  °C 
per  kilometer  of  height  above  the  earth's  surface.  The  energy  flow  rate  in 
the  convection  column  above  a  line  of  fire  (Pf)  is  the  rate  at  which  thermal 
energy  is  converted  to  kinetic  energy  at  the  height  specified  for  Pw  in  the 
convection  column.  Byram  pointed  out  that  his  case  studies  of  several  wild- 
fires had  shown  that  the  ratio  Pf/Pw  can  be  useful  in  understanding  and 
predicting  the  onset  of  erratic  fire  behavior  and  the  occurrence  of  blowup 
fires.  According  to  Byram,  his  studies  demonstrated  a  strong  relationship 
between  the  occurrence  of  blowup  fires  and  values  of  Pf/Pw  equal  to  or 
greater  than  1  for  at  least  1,000  feet  above  the  fire. 

Other  fire  researchers  have  used  the  Pf/Pw  concept.  It  has  been  used  in 
Australian  and  Russian  studies  of  spotting  and  crown  fire  spread  (Cheney 
and  Bary  1969;  Kurbatsky  1969;  Mateev  and  Mateev  1977).  In  the  United 
States,  Wade  and  Ward  (1973)  and  Simard  and  others  (1983)  used  Pf/Pw 
in  the  analysis  and  interpretation  of  wildfire  behavior.  Aronovitch  (1989) 
suggested  using  PflPw  to  forecast  the  occurrence  of  blowup  and  to  estimate 
the  height  of  the  convection  column.  He  recommended  that  wind  profiles 
needed  for  computing  Pw  be  made  a  part  of  fire  management  systems  such 
as  BEHAVE.  Rothermel  (1991a)  developed  a  chart  for  estimating  the  rate 
of  spread,  rate  of  energy  release,  and  flame  length  of  crown  fires  in  the  North- 
ern Rocky  Mountains.  He  then  used  Pf/Pw  to  predict  whether  a  given  fire 
will  be  wind-driven  (with  a  plume  tipped  away  from  vertical  in  the  direc- 
tion of  the  wind)  or  plume-dominated  (with  a  convection  column  standing 
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nearly  vertical).  In  addition,  Byram  and  others  (1964)  compared  Pf/Pw  with 
observations  of  flame  characteristics  in  a  laboratory  wind  tunnel  study.  More 
recently,  Martin  and  others  (1991)  used  Pf/Pw  to  study  how  flame  angle  is 
related  to  windspeed  for  laboratory  pool  fires,  fires  in  which  liquid  fuels  are 
burned  in  pans  or  trays. 

With  the  increasing  interest  in  using  Pf/Pw  to  interpret  and  characterize 
the  behavior  of  wildfires  and  smaller  scale  fires,  it  is  likely  that  fire  manag- 
ers and  fire  researchers  will  become  aware  of  an  inconsistency  in  the  units 
of  Pw  and     in  Byram's  original  work  (1959)  and  in  later  documentation 
(Brown  and  Davis  1973).  The  problem  is  as  follows:  units  of  Pw  and  Pf  are 
given  as  ft  lb/ft2  s,  where  lb  must  refer  to  pounds  force  to  obtain  a  flow  of 
energy.  Substitution  of  units  into  the  equations,  however,  shows  that  lb  in 
the  numerator  of  Pw  and  Pf  has  units  of  pounds  mass  rather  than  pounds 
force.  This  inconsistency  may  be  a  source  of  confusion  for  persons  not  famil- 
iar with  the  systems  of  units  commonly  employed  in  the  physical  sciences. 
The  problem  also  becomes  evident  when  one  wishes  to  calculate  Pw  and  Pf 
separately  in  units  differing  from  those  given  by  Byram. 

This  note  (1)  demonstrates  the  inconsistency  in  units  of  the  original  Pw 
and  Pf  equations,  (2)  introduces  a  general  version  of  the  equations  valid  in 
any  system  of  units,  (3)  describes  conversion  factors  required  for  calculating 
Pw  and  Pf  in  the  British  gravitational  and  SI  systems,  and  (4)  presents  sam- 
ple calculations  in  the  British  and  SI  systems  using  the  general  equations. 

ORIGINAL  EQUATIONS  FOR  Pw  AND  Pf 

The  original  equations  for  Pw  and  Pf  in  a  neutrally  stable  atmosphere  as 
given  by  Byram  (1959)  are 


where  Pw  is  the  rate  of  energy  flow  in  the  wind  field  at  height  z  above  a  fire 
and  Pf  is  the  rate  at  which  thermal  energy  is  converted  into  kinetic  energy 
in  the  convection  column  at  height  z.  The  quantities  in  equations  (1)  and 


(2)  and  their  units  in  the  British  system  are 

P 

=  air  density  at  height  z,  lbm/ft3 

V 

=  windspeed  at  height  z,  ft/s 

r 

=  forward  rate  of  spread,  ft/s 

g 

=  acceleration  of  gravity,  ft/s2 

I 

=  fireline  intensity,  Btu/ft-s 

cp 

=  specific  heat  of  air  at  constant  pressure,  Btu/lbm-°F 

T 

x  0 

=  free-air  temperature  at  the  elevation  of  the  fire,  °F. 

In  the  strictest  sense,  units  of  energy  in  the  British  system  are  ft  lbf,  but  the 
energy  units  divide  out  cleanly  in  the  ratio  I/cp  of  equation  (2).  Thus,  the 
more  familiar  Btu  of  energy  is  used  here.  Substitution  of  units  into  equations 
(1)  and  (2)  shows  that  both  Pw  and  P^have  units  of  lbm/ft-s,  or  ft  lbm/ft2-s 
when  multiplied  by  the  ratio  ft/ft.  Only  for  a  system  of  units  in  which  1  lbf = 
1  lbm  will  equations  (1)  and  (2)  provide  numerically  correct  answers.  In 
general,  the  equations  are  not  formulated  with  the  anticipation  of  possible 
use  in  other  systems  of  units. 

The  inconsistency  in  units  is  related  to  Newton's  second  law  of  motion, 
F  =  ma  =  mg,  and  to  the  use  of  absolute  or  gravitational  units  in  the 


Pw  =  p(v  -  rf/2g 
Pf  =  I/cp(T0  +  459) 


(1) 


(2) 
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equations  for  Pw  and  Pf.  In  Newton's  law,  F  denotes  a  force,  m  is  a  mass,  a  is 
an  acceleration,  and  g  is  the  acceleration  due  to  gravity.  In  an  absolute  system, 
such  as  the  SI  system  of  units,  there  is  one  unit  for  each  physical  quantity 
and  the  product  mg  equals  the  force  F  by  definition  of  the  units.  In  the  gravita- 
tional system,  however,  mass  is  divided  by  a  conversion  factor  that  is  nu- 
merically equal  tog  but  has  different  units.  This  factor  is  denoted  here  by 
g0,  and  Newton's  law  becomes  F  =  mglgQ.  Thus,  a  force  of  1  pound  (1  lbf)  is 
numerically  equivalent  to  a  mass  of  1  pound  (1  lbm).  When  g  has  its  stan- 
dard value  of  32.2  ft/s2,  multiplication  by  the  ratio  g/g0  renders  each  equa- 
tion consistent  with  respect  to  units  because  the  conversion  factor  is  given 
by  g0  =  32.2  ft  lbm/lbf  s2.  Thus,  g/g0  numerically  equals  1,  and  equations 
(1)  and  (2)  are  better  expressed  as 

Pw  =  p(v  -  r)V2g0  (3) 

Pf  =gI/g0cp(T0  +  459).  (4) 

Substitution  of  the  appropriate  units  now  leads  to  ft-lbf/ft2-s,  the  units  for 
Pw  and  Pf  given  by  Byram  (1959).  For  consistent  units,  equations  (3)  and 
(4)  should  be  used  in  place  of  equations  (1)  and  (2)  when  calculations  of  Pw 
and  Pf  are  made  in  the  British  gravitational  system. 

GENERAL  EQUATIONS  FOR  Pw  AND  Pf 

Equations  for  Pw  and  Pf  may  be  written  in  a  form  applicable  to  any  system 
of  units.  Transformation  from  one  system  to  another  is  easily  accomplished 
if  a  table  of  conversion  factors  relating  mass,  force,  and  energy  is  available. 
The  general  equations  correspond  to  equations  (3)  and  (4),  with  g0  dimen- 
sionless  and  equal  to  1.  In  addition,  an  absolute  temperature  is  required  in 
place  of  (T0  +  459)  in  equation  (4).  The  equations  are 

Pw  =  p(v  -  rf/2  (5) 

Pf  =gI/cpTa  (6) 

where  Ta  is  the  absolute  temperature  of  air  at  the  elevation  of  the  fire  in  the 
units  employed  and  all  other  quantities  have  been  defined  previously.  Equa- 
tions (5)  and  (6)  correspond  to  the  Pw  and  Pf  relationships  derived  by  Byram 
from  consideration  of  the  physical  processes.  I  am  writing  an  article  pre- 
senting the  derivation  of  Pw  and  Pf. 

Factors  are  needed  for  converting  from  force  to  mass  times  acceleration 
and  from  force  to  energy  per  unit  length.  For  the  British  and  SI  systems 
of  units  these  factors  are 

British:  1  lbf  =  32.2  lbm-ft/s2  =  32.2  lbm-ft2/ft-s2  (7) 

SI:  1  newton  =  1  kgm/s2  =  U/m  =  0.001  kWs/m.  (8) 

These  identities,  when  used  with  equations  (5)  and  (6),  facilitate  calculation 
of  Pw  and     in  units  of  energy  per  unit  time  per  unit  area,  or  equivalently, 
power  per  unit  area  (Martin  and  others  1991;  Rothermel  1991b). 

SAMPLE  CALCULATIONS 

Suppose  that  at  some  specified  height  in  the  atmosphere  above  a  line  fire  the 
quantities  p  and  v  have  values  in  British  units  of  0.0695  lbm/ft3  and  20  ftVs, 
respectively.  Further  suppose  the  forward  rate  of  fire  spread  is  2  ftVs  and  the 
corresponding  fireline  intensity  is  15,000  Btu/ft-s.  At  the  elevation  of  the  fire, 
cp  and  Ta  have  values  of  0.24  Btu/lbm-°R  and  540  °R.  These  quantities  are 
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used  with  tables  of  conversion  factors  to  compute  Pw,  Pf,  and  the  ratio  Pf/Pw 
in  the  British  and  SI  systems;  the  calculations  are  given  below.  The  SI  cal- 
culations are  checked  by  conversion  back  to  the  British  system,  illustrating 
the  transformation  from  one  system  of  units  to  another. 

British:  From  equations  (5)  and  (7), 

Pw  =  [0.0695  lbm/ft3  ■  (20  -  2)3  ft3/s3][l  lbf.ft-s2/32.2  lbm-ft2]/2  =  6.30  ft-lbf/ft2*. 

From  equations  (6)  and  (7), 

Pf  =  [32.2  ft/s2  •  15,000  Btu/ft-s][l  lbf-ft-s2/32.2  lbm-ft2]/[0.24  Btu/lbm-0R  •  540  °R] 
=  115.74  ft-lbf/ft2-s 

and  Pf/Pw  =  115.74/6.30  =  18.37. 

SI:  Suppose  the  fire  and  atmospheric  variables  presented  earlier  are  avail- 
able in  SI  units.  Then 


p 

=  1.11  kg/m3 

CP 

=  1,003.9  J/kg-K 

Ta 

=  300  K 

V 

=  6.1  m/s 

r 

=  0.61  m/s 

I 

=  5.19-  104kW/m 

g 

=  9.8  m/s2. 

From  equations  (5)  and  (8), 

Pw  =  [1.11  kg/m3  •  (6.1  -  0.61)3  m3/s3][0.001  kW-s/m/1  kg-m/s2]/2  =  0.0918  kW/m2. 

From  equations  (6)  and  (8), 

Pf  =  [9.8  m/s2  ■  5.19  •  104  kW/m][l  J/m/1  kg-m/s2]/[l,003.9  J/kg-K  •  300  K]  = 
1.69  kW/m2 

and  Pf/Pw  =  1.69/0.0918  =  18.41,  in  close  agreement  with  the  computation 
using  the  British  system. 

Values  of  Pw  and     in  SI  units  are  converted  to  British  units  as  follows: 

Pw  =  0.0918  kW/m2  ■  1  ftlbf/0.001356  kWs  •  1  m2/10.76  ft2  =  6.29  ft  lbf/ft2-s 

where  use  is  made  of  the  identities  1  ft-lbf  =  1.356  J  =  0.001356  kW-s  and 
1  m2  =  10.76  ft2. 

Similarly, 

Pf  =  1.69  kW/m2  •  1  ft-lbf/0.001356  kWs  •  1  m2/10.76  ft2  =  115.83  ft  lbf/ft2-s. 

The  negligible  difference  in  Pw  and  Pf  values  between  the  two  systems  is 
due  to  rounding  errors. 
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